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Electron-electron interactions are responsible for a correction to the conductance of a diffusive 
metal, the "Altshuler-Aronov correction" 6G aa- Here we study the counterpart of this correction for 
a ballistic conductor, in which the electron motion is governed by chaotic classical dynamics. In the 
ballistic conductance, the Ehrenfest time te enters as an additional time scale that determines the 
magnitude of quantum interference effects. The Ehrenfest time effectively poses a short- time thresh- 
old for the trajectories contributing to the interaction correction. As a consequence, SG aa becomes 
exponentially suppressed if the Ehrenfest time is larger than the dwell time or the inverse tem- 
perature. We discuss the explicit dependence on Ehrenfest time in quasi-one and two-dimensional 
antidot arrays. For strong interactions, the sign of SGaa may change as a function of temperature 
for temperatures in the vicinity of H/tb- 
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I. INTRODUCTION 

Electronic transport in weakly disordered metals 
is successfully described by the Boltzmann theory, 
in which electrons are treated as effectively classi- 
cal particles moving freely between scattering events. 
The wave nature of electrons gives rise to a num- 
ber of corrections to transport properties, such as the 
weak localization correction, ^'^ the Altshuler-Aronov in- 
teraction correction, ^'^ or the universal conductance 
fluctuations.^'^ Weak localization results from the con- 
structive interference of electrons propagating along 
time-reversed paths. ^ The physical intuition behind the 
interaction correction is constructive interference of elec- 
tron trajectories which are scattered on impurities and 
Friedel oscillations of the electron density. ^'^ These quan- 
tum corrections become increasingly important as the 
temperature is lowered, the effective dimensionality of 
the sample is reduced, or as the disorder level is increased. 
They have a distinctive and universal dependence on ex- 
ternal parameters, such as temperature or magnetic field, 
which makes them identifiable in experiments. In partic- 
ular, the two quantum corrections to the conductivity, 
weak localization and the Altshuler-Aronov correction, 
can be distinguished by application of a magnetic field, 
since weak localization is suppressed by already a very 
small magnetic field, whereas the Altshuler-Aronov cor- 
rection is not. 

A "classical analog" of a disordered metal is realized 
in high-mobility semiconductor structures with randomly 
placed large antidots.^^'^^ The absence of impurities en- 
sures that electrons move ballistically between reflections 
off the antidots. The reason why these systems are re- 
ferred to as classical is that the size of the antidots a is 
much larger than the Fermi wavelength Ap. As a re- 
sult, not only the electron's motion through the two- 
dimensional electron gas, but also the reflection off an 
antidot is described by classical mechanics. (In contrast, 
in a disordered metal, the size of impurities is compara- 



ble to Ap, so that the scattering event is strongly diffrac- 
tive.) For an irregular arrangement of antidots, the clas- 
sical dynamics is chaotic. Nearby trajectories separate 
exponentially in time, the exponential separation being 
characterized by the Lyapunov coefficient A. The chaotic 
dynamics is essential for the existence of quantum cor- 
rections in this system, as it magnifies the quantum un- 
certainty of even a minimal wavepacket up to classical 
dimensions within the short time 

Te = ^ln(a/AF), (1) 

thus transforming the classical dynamics into quantum- 
diffract ive dynamics on time scales larger than te-^^ The 
time Te is known as the "Ehrenfest time" . 

Since wave effects are not operative for times shorter 
than Te — electrons essentially move along classical tra- 
jectories up to the Ehrenfest time — , the Ehrenfest time 
serves as a short-time threshold for the duration of the 
trajectories contributing to the quantum corrections in 
an antidot array. For weak localization, it was found 
that the correction to the conductivity is exponentially 
suppressed if te is larger than the dwell time td , the typ- 
ical time to be transmitted through the system, or the 
dephasing time r^.^^^^ In contrast, other quantum cor- 
rections, such as the universal conductance fluctuations, 
remain finite if te ^ td."*^^'^^ ""^^ 

The goal of this article is to present a theory of 
the Ehrenfest-time-dependence of the Altshuler-Aronov 
correction SGaa- Our analysis significantly extends a 
previous calculation by Kupferschmidt and one of the 
authors, which studied the te dependence of the inter- 
action correction to the conductance of a ballistic double 
quantum dot and found that 6G aa is strongly suppressed 
if Te exceeds the dwell time td or the inverse tempera- 
ture h/T. The double quantum dot studied in Ref. 19 
is the simplest system with nonzero Altshuler-Aronov 
correction to the conductance, and is characterized by 
a long-range interaction, which is spatially homogeneous 
within each dot. The theory presented here is valid for 
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both short-range and long-range interactions and can be 
apphed to any geometry in which the classical electron 
dynamics is chaotic — although we will focus our dis- 
cussion on the case of an antidot array. For the general 
case considered here we confirm the suppression of SG aa 
for te ^ min(rD, ^/T) and we calculate the precise func- 
tional dependence of SG aa on td and T for finite Ehren- 
fest time. The explicit dependence on temperature is 
characteristic for the interaction correction, which has its 
origin in virtual processes with an energy transfer larger 
than temperature. 

Our calculation makes use of a semiclassical formalism 
that starts from the saddle-point approximation around 
classical trajectories for the single-particle Green func- 
tion. In this way, the conductance in the absence of 
electron-electron interactions is written as a double sum 
over classical trajectories that connect source and drain 
reservoirs. ^^'^^ Weak localization and other quantum cor- 
rections to the conductance then follow from special con- 
figurations of trajectories, in which the two trajectories in 
the summation are piecewise paired, and proceed through 
"crossings" at points where the pairing is changed. ^^'^^ In 
the language of diagrammatic perturbation theory, seg- 
ments where the trajectories are paired correspond to dif- 
fusons or cooperons, whereas the crossings correspond to 
Hikami boxes. The application to interacting electrons 
requires a modification of the formalism, which will be 
described in detail below. 

Our analysis applies to a "ballistic" conductor, where 
the label "ballistic" is meant to specify that the elec- 
trons move along well-defined classical trajectories. In 
the literature, "ballistic" sometimes refers to a different 
limit, and several calculations of the interaction correc- 
tion to the conductance have been reported for such "bal- 
listic limits" . Whereas the original work of Altshuler and 
Aronov^ addressed a disordered metal with short-range 
scatterers in the diffusive regime Tr <C 1, the theory was 
generalized to account for the effects of higher temper- 
atures Tr ^ 1, a regime referred to as "ballistic" .^'^^'^^ 
The case of a smooth disorder potential, in which scat- 
tering is predominantly forward, was considered in Ref. 
26. Another type of system, where interaction correc- 
tions appear, are networks of capacitively coupled ballis- 
tic quantum dots,^^^^ where, however, Ehrenfest-time- 
related phenomena can be neglected as long as te is much 
smaller than the dwell time in a single quantum dot. 

In Section II we present our theory of the Ehrenfest- 
time dependence on the interaction correction for a 
generic ballistic chaotic conductor. In Sec. Ill we then ap- 
ply our formalism to an antidot array, where the classical 
electronic motion is diffusive on length scales much larger 
than the spacing between antidots, and the Coulomb in- 
teraction is dynamically screened by the diffusively mov- 
ing electrons. For the antidot array, we find 6G aa oc 
exp(— te/td — ^-kTt^/K) in the limit that the Ehrenfest 
time is larger than dwell- time and inverse temperature. 
For small Ehrenfest times, we recover the results of a 
disordered metal with quantum impurities, which show a 



much weaker temperature dependence (algebraic or log- 
arithmic, depending on dimensionality). We conclude in 
IV. 



II. SEMICLASSICAL THEORY OF THE 
INTERACTION CORRECTION 

In this section we present the semiclassical descrip- 
tion of the interaction corrections for a conductor with 
a well-defined chaotic classical electron dynamics. We 
first review the expressions for the interaction corrections 
to the conductance in terms of the single-particle Green 
function, and then apply the semiclassical approximation 
methods, taking into account the finite Ehrenfest time. 



A. Skeleton diagrams for the conductance 

For definiteness, we consider a two-dimensional ballis- 
tic conductor, such as a ballistic electron gas with an 
antidot array, in contact with reservoirs at x = and 
X = L, see Fig. 1. Without interactions, we can calculate 
the conductance G from the Kubo formula, 

X [v,g^{v,v'-i)v,,g^(y' ,v-0].'=o , (2) 

x = L 

where /(^) = l/(exp(^/T) + 1) denotes the Fermi func- 
tion, 

is the velocity operator, and S^(r, r'; ^) and S^(r, r'; ^) is 
the retarded and advanced single-particle Green function, 
respectively. Retarded and advanced Green functions are 
related as 

gA(r',r;0 = a^(r,r';0*- (4) 

To leading (first) order in the interaction strength, 
the interaction correction 5Gaa is obtained by replac- 
ing g\v, r'; by C?^(r, r'; i) + 5gf{v, r'; O + '^^h (r, r'; C) 
and expanding to first order in the interaction [/^9'i9,3i 
where 

5g^{v,v',^) = J^J dv.dr, tanh(^^) 

xG^(r,ri;0G'^(r2,r';0 
X {C/A(ri,r2;w)G^(ri,r2;C-w) 
-C/i^(ri,r2;:^)GA(ri,r2;^-a;)}, 

(5) 
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FIG. 1: (Color online) Schematic picture of the system under 
consideration: A ballistic conductor, attached to ideal leads 
at X = and x — L. In the semiclassical calculation of the 
conductance, one retarded and one advanced Green function 
are attached to a current vertex located at the interface with 
the leads. In a semiclassical picture, these Green functions 
are associated with "retarded" and "advanced" classical tra- 
jectories (solid and dashed in the figure), both of which must 
point into the conductor (1). A current vertex combined with 
two Green functions of the same kind is not possible for the 
calculation of the conductance: after pairing, we have trajec- 
tories that go straight into the leads (2). On the contrary, for 
a calculation of the conductivity, the current vertex can be 
anywhere inside the conductor, and pairing of retarded and 
advanced trajectories is possible also if two Green functions 
of the same kind are attached to one current vertex (3 and 
4). 

xG^(r,ri;OG^(ri,r';C) 
X {C/A(ri,r2;0)G^(r2,r2;e-t^) 
-f/f^(ri,r2;0)GA(r2,r2;^-w)}, 

(6) 

and with similar expressions for the advanced func- 
tions SQ^{r\ r; ^) and SQ^{r\ r; In these expressions, 
[/^(ri, r2; cj) and /7^(ri, r2; cj) are the retarded and ad- 
vanced interaction kernels, respectively. The interaction 
is taken to be zero in the leads, for x < and x > L. Such 
a structure represents the change of the single-particle 
Green function due to scattering off Friedel oscillations 
of the density matrix p(ri, (Fock) or the density p{t2) 
(Hartree).^ The resulting contributions to SG aa are rep- 
resented diagrammatically as in Fig. 2. 

We would like to emphasize that we kept here only 
those diagrams for which one retarded and one advanced 
Green function are connected to each current vertex. 
These are the relevant diagrams for the calculation of 
the conductance. This is in contrast to the calculation 
of the interaction correction to the conductivity 5a aa^ 
where diagrams with two Green functions of the same 
kind attached to the current vertex play an important 
role.^ (The conductivity a is expressed in a similar way 
as Eq. (2), but contains integrals over the x-coordinates 
as well, rather than fixing them to the contacts at x = 




FIG. 2: (Color online) Hartree (upper) and Fock (lower) di- 
agrams: solid (dashed) lines represent retarded (advanced) 
Green functions, the wiggly line represents the interaction. 
Each diagram has a counterpart with retarded and advanced 
Green functions interchanged. 



and X = L). The difference between the conductance cal- 
culation and the conductivity calculation is readily seen 
in the semiclassical language. In that language. Green 
functions are associated with classical trajectories, and 
only terms in which "retarded" and "advanced" trajecto- 
ries are paired contribute. (For more details, see below). 
Since the leads are assumed to be free of disorder and 
without electron-electron interactions, both the retarded 
and the advanced trajectories at the positions r and in 
Eq. (2) must point into the conductor if the conductance 
is calculated. On the other hand, for a current vertex in 
the system's interior, pairing of advanced and retarded 
trajectories is still possible even if two Green functions 
of the same kind are attached to the current vertex, see 
Fig. 1. 

The fact that different diagrams are needed for the cal- 
culation of conductance and conductivity is well known, 
the same is true for the Drude conductance and conduc- 
tance fluctuations of a disordered metal (see Ref. 32). For 
instance, for the calculation of the Drude conductivity ctq 
in a metal with short-ranged disorder there is no need to 
dress the diagram with an impurity ladder, while for the 
classical conductance Go, the diagram dressed with an 
impurity ladder, i.e., a diffuson, is most relevant. For the 
Drude conductivity one might argue that the distance be- 
tween the current vertices is of the order of the mean free 
path, since the two Green functions decay on this scale. 
The diffuson in turn is long-ranged, and hence needed 
to describe propagation from one lead to the other, as 
required for the conductance. 

We also note that the calculation of the conductance 
as it is outlined here is similar to the calculation of the 
density-density correlation functions. ^^'^^ A subtle point 
in this regard is the existence of additional corrections 
in the calculation of the density-density correlation func- 
tion, namely vertex corrections and the so-called wave- 
function renormalization. In the description developed 
below, both of them appear to vanish. For the density- 
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density correlation function, in turn, vertex corrections 
and the wave function renormalization cancel each other, 
so that these corrections do not lead to a net change of 
the result in either case. 



B. Semiclassical theory 

The conductance G depends on the precise locations of 
antidots, the system boundary, and on the Fermi energy. 
We now employ a semiclassical analysis in order to iden- 
tify those contributions to the conductance that remain 
after an average over the Fermi energy. 

Starting point is the semiclassical expression of the 
Green function ^^(r, r';^) as a sum over classical tra- 
jectories a from to r at energy 



2-K 



(7) 



Here, <Sa(r, r';^) is the classical action corresponding to 
the trajectory a, which has the properties 



dv 



Pa, 



dr' 



Pa. 



and 



(8) 



(9) 



where and denotes the momentum at the end and 
beginning of ce, respectively, and Tq, is the duration of 
the traj ectory. Th e stability amplitude Aq, is given by 
Ao, = ^/\ det{Da)\, with 



dr'dr 



(10) 



The semiclassical Green function further contains an ad- 
ditional phase-shift, the so-called Maslov index, which 
we omitted because it does not play a role in our con- 
siderations. The semiclassical expression for advanced 
Green function follows from Eq. (4). 

Using the semiclassical Green function (7) we express 
the interaction correction SG aa as a fourfold sum over 
classical trajectories. We refer to these trajectories as 
"retarded" or "advanced" , depending on the type of the 
Green function that they originate from. The summation 
over classical trajectories can be simplified for a system 
with chaotic classical dynamics: In this case, the classi- 
cal trajectory and hence the classical action depend very 
sensitively on the initial conditions. On the other hand, 
in the semiclassical limit h ^ only configurations of 
trajectories with sum of the actions of the "retarded" tra- 
jectories systematically equal to the sum of the actions of 
the "advanced" trajectories up to a difference AS of the 
order of h contribute substantially to the conductance. 
This occurs only if the "retarded" and "advanced" tra- 
jectories are piecewise paired, whereby they can exchange 



"partners" only at a "small-angle encounter" at which 
two pairs meet to within a phase-space distance of order 

For the remaining summation over trajectories, we use 
a sum rule that expresses the summation over trajectories 
a between positions and r and at energy ^ in terms 
of an integral over the trajectory's duration t, the initial 
and final momenta and p, as well as a "trajectory 
density" p^(X' X; t) between the phase-space points 
X' = (r^pO and X = (r,p), 36,37 



(11) 



see App. A 1 for details. Here / is an arbitrary func- 
tion. The initial point in phase space X' = (r',p'), 
together with the Hamilton function H uniquely deter- 
mines the classical trajectory, and after a time t this 
trajectory has reached the phase space point X(t) = 
(r(r', p'; t), p(r', p', t)). The trajectory density 

/)(X'^X;t) = J[X-X(t)] (12) 
= ^[r-r(r^p^^)Hp-p(r^p^^)] 

selects then only those phase space points which are con- 
nected by a trajectory of duration t. At fixed energy 
^, the momentum integrations are restricted to the en- 
ergy shell, dp^ = dpS{^ — H{p, r)), and for the trajectory 
density we may factor out the part which ensures energy 
conservation, 

p{X'^X;t) = p^{X' ^X;t)S[H{X)-H{X')].{13) 

The factor provides the Jacobian for this transforma- 
tion. 

Following the procedure outlined so far, we obtain an 
expression in terms of trajectory densities which, strictly 
speaking, is a sum of ^-functions. We then replace the 
exact trajectory density by a coarse-grained smooth 
density p.^^'^^ The coarse graining takes place with re- 
spect to small fiuctuations of the initial and final phase 
space points and/or the positions of the scattering discs 
or the system's boundaries. In the regime Atd ^ 1, 
where the chaotic dynamics has fully developed, the clas- 
sical dynamics is essentially stochastic, which justifies the 
coarse graining procedure. The coarse-grained trajectory 
density 



p^(X'^X;i) = P(X,X';t) 



(14) 



can be identified with the probability density P(X, X'; t) 
that a particle originating at the phase space point X' = 
(r', p) is found at the phase space point X = (r, p) after 
a time t. (Since we are interested in the regime where 
temperature is much smaller than Fermi energy, we drop 
the dependence of the classical propagators on ^.) For 
the case of antidot arrays, this probability density is de- 
scribed by a diffusion equation. 
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We now turn to the semiclassical calculation for the in- 
teraction correction. 



FIG. 3: (Color online) Configuration of trajectories relevant 
for the first diagram to the Fock contribution. 



The Drude conductance is obtained by keeping only 
pairs of classical trajectories that connect source and 
drain reservoirs. Following the steps described above, 
we find 



[^,P(X,XV)<]. 



(15) 



Fock contribution 



We start with the Fock contribution to SG aa, which is 
given by the two lower diagrams in Fig. 2, together with 
their counterparts, which are obtained by interchanging 
the retarded and advanced Green functions. We first 
consider the conductance correction SG^)^ from the lower 
left diagram and its counterpart, which reads 



X [i),g^(r, r'; ^^.'^^(r', ra; Oa^(r2, ri; ^ - uj)g^{ri,r; 0] x'=o | • (16) 

x = L I 

After insertion of the semiclassical expression Eq. (7) for the Green functions, we obtain a sum over one retarded and 
three advanced trajectories. In the semiclassical limit, a convolution of Green functions is customarily calculated using 
the stationary phase approximation. For this one first needs to determine the configurations of trajectories which 
make the total action stationary. This results in a factor ex 6**^^% where Sst is obtained by inserting the stationary 
configuration into the total action. Integration over quadratic fluctuations around the stationary configurations then 
renders the prefactor. In the present case the convolution of Green functions is accompanied by the interaction 
propagator. One might expect, that the interaction propagator affects the stationary trajectories, such that they no 
longer connect to a single classical trajectory, as in the case without interaction. However for the calculation of the 
conductance, we need to pair the advanced trajectories with the retarded one, see Fig. 3. Hence, performing the 
integration over ri and r2 in Eq. (16) within stationary phase approximation, we only take into account stationary 
configurations that connect to a single classical trajectory. The detailed calculation is carried out in Appendix A 2 
and has the result 
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V A^e-^-^^/ft r°dt /'*d^'[/R(r«(^),r„(^');c^)e^'^(*-*')/^ (17) 

-^0 Jo 



t9 {-2mhf/^ 

where rQ,(t) is the coordinate of trajectory a after time t. The integration over time reflects the freedom to choose 
ri and anywhere along the trajectory a; the factor e^^^^~^ takes into account the action difference at different 
energies, ^^(^ — cj) = 5^(0 ~ ^^a. (foi" cj <C 0- With XQ,(t) = (rQ,(t), PQ;(t)) we may rewrite 

[/^(r,(t),r,(tO;a;) = ^ dXidX2p^(X:, ^ Xi; t)p^(X:, ^ X2; 0/7^(ri, r2; c^), (18) 

where (iX = drdp^ is an integration over phase-space points on the energy shell, and X'^, = Xq;(0) is the initial 
phase-space point of trajectory a. 

After inserting Eqs. (17) and (18) into Eq. (16), and upon applying the semiclassical approximation to the retarded 

F 1 

Green function as well, the interaction correction SGj^j^ is expressed as a double sum over trajectories a and P running 
from to r. Only diagonal terms with a = ^ are systematically nonzero, so that we only keep these. Again making 
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use of the sum rule (Eq. (11)) we express ^G^'a integral over the two intermediate phase space points Xi 

and X2. Before the exact trajectory densities can be replaced by their coarse-grained versions, we split the classical 
trajectories into uncorrelated segments using the equality 



Jo Jo Jo Jo 



(19) 



After coarse-graining, the expression for ^G^'a involves the probability densities P(Xi X2;t) for the chaotic 
classical motion. The expression can be further simplified by introducing 

Pi„(X) = jdy' J dp'^ dt KP(r', p' ^ X; t)],,^^ , 



Pout (X) = J dy j dp^ J dt [P(X r, p; i)ux]^=L , 



(20) 



which express the probability that a trajectory at phase space point X entered at the left contact or exits at the right 
contact, respectively. Using the equality 



we finally obtain 



SGi'l = 



where we singled out the part containing classical propagators 



J dw-^ (wcoth^)lm|y" dXiC?X2[/^(ri,r2;w)i^i(Xi,X2;w)|, (22) 



POO 

ifi(Xi,X2;w)= / dtPout(Xi)P(Xi,X2;i)e^"*/'^Pi„(X2). 

^0 



(23) 



We now consider the interaction correction SG^^ that corresponds to the lower right diagram of Fig. 2 and its 
counterpart obtained by switching retarded and advanced labels, 

= - v/'^^ (-^) / S*""^ (^) '"^l/'^^Z ""y'J dridr2U^ir,,r2,.) 

X [t).a^(r, n; C)a^(ri, r2; ^ - o:)g^{r2, r'; 0«.'^^(r', r; 0] | • (24) 

x = L I 



Insertion of the semiclassical expression for the Green 
functions leads to a fourfold sum over two retarded tra- 
jectories (from to r2 and from ri to r), and two ad- 
vanced trajectories (from to r and from ri to r2). 

Because of the specific requirements for the start and 
end points of the trajectories, it is not possible, to pair 
the trajectories one by one for their entire duration. In- 
stead, the trajectories need to undergo a "small-angle 
encounter", in which all four trajectories are close to- 
gether in phase space for at least part of their length. 
The four possible configurations of trajectories are shown 
in Fig. 4, where we take into account the possibilities 
that none, one, or both points ri and r2 lie inside the 
encounter region. Their contributions to SGaa will be 
denoted SG^^l^-SG^^l^ , see Fig. 4. 

The summation over classical trajectories with a small- 
angle encounter follows the procedure outlined in Refs. 



14,23. We refer the reader to appendix A3 for details, 
and proceed with the results of that calculation. All four 
contributions to SGaa have the same form as the con- 
tribution from the first diagram, see Eq. (22), but with 
different expressions for the function i^(Xi, X2; cj). For 
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FIG. 4: (Color online) Configurations of trajectories relevant 
for the second diagram to the Fock contribution. Encounter 
regions are indicated in blue; here the motion of all four tra- 
jectories is correlated. 



the contributions 2a-2d these expressions read 



-J dXcZX'Pi„(X)P(X,Xi;w) 

xP(X2,X';w)Pout(X') 
d 



9tf 



(25) 



(26) 



P(X',X;TE)e'"^/'' 

-y"dXPi„(X)P(X,Xi;w)Po„t(X2) 
xP(X2,X;TE)e'""^/^ 

- J dX'Pin(Xi)P(X2,X»Pout(X') 

X P(X',Xl;TE)e*'^^/^ (27) 

- PdiPin(Xi)P(X2,Xi;t) 
Jo 

X e*"*/''Pout(X2). (28) 
Here P(X, X'; w) is the Fourier transform of P(X, X'; t), 



P(X,X';a;)= / dtP(X, X'; t)e*"*/''. 

JO 



(29) 



Taken together, Eqs. (22), (23), (25), (26), (27), and (28) 
determine the general result for the Fock contribution 
to SGaa for finite Ehrenfest time, expressed in terms of 
classical propagators. 

Let us briefly discuss the effect of the Ehrenfest time: 
Interestingly, the contribution does not involve a 
crossing and therefore shows no dependence on the 
Ehrenfest time. However, it is cancelled by the contri- 
bution i^^^, if the travel time between Xi and X2 is 
shorter than the Ehrenfest time. Thus, adding all con- 
tributions together, we indeed find, that effectively only 
trajectories with a duration longer than the Ehrenfest 
time are responsible for the interaction correction. 



D. Hartree contribution 

The Hartree contribution to the Altshuler-Aronov cor- 
rection is given by the two upper diagrams in Fig. 2. 
Proceeding as in the case of the Fock contribution, each 
Green function is written as a sum over classical trajec- 
tories, which must then be piecewise paired in order to 
give a nonvanishing contribution to the interaction cor- 
rection to the conductance. The resulting configurations 
of classical trajectories are shown schematically in Fig. 
5. The trajectory configurations of Fig. 5 are in one- 
to-one correspondence to those of Figs. 3 and 4 for the 
Fock contribution to SGaa- The diagram of Fig. 5a cor- 
responds to that of Fig. 3, whereas the diagrams of Fig. 
5b-e correspond to those of Fig. 4a-d. 

Unlike the Fock diagrams, all diagrams for the Hartree 
correction involve a finite-angle crossing of the trajecto- 
ries, in addition to the small-angle encounter of Figs. 
5b-e. Another important difference is that the action 
difference AS for the Hartree case depends on the two 
positions ri and r2 associated with the interaction vertex. 
Denoting the momenta involved in the finite-angle cross- 
ing of the trajectories by pi and p2, see Fig. 5, the action 
difference contributes an additional oscillating phase fac- 
tor e*(^i-^2) (pi-P2)/^ (No fast oscillating phase factors 
are associated with the integration over ri and T2 for the 
Fock diagrams.) For chaotic classical motion, the direc- 
tions of the momenta pi and p2 are random and uncor- 
related, while the magnitude |pi| = |p2| = Pf is fixed 
by energy conservation. As a result, only a short-range 
component of the interaction contributes to the Hartree 
correction, and one finds the same expression for SGaa 
as for the Fock contribution, with the replacement^ 

[/^(ri-r2;a;) ^ -2(5(ri - r2) (30) 

where the brackets (...) indicate an average over the mo- 
menta pi and p2 with |pi| = |p2| = pp. In case of a 
short-range interaction /7(ri,r2) cx S{ri — r2), one veri- 
fies that this replacement rule leads to SG^j^ = — 2(5G^A5 
in agreement with Eqs. (5) and (6). 

Intuitively, the interaction correction associated with 
the trajectory configurations of Fig. 5 can be inter- 
preted as the interference of electrons that follow a 
classical trajectory connecting source and drain con- 
tacts, and electrons that additionally scatter from Friedel 
oscillations. ^'^'^^ In the configuration of Fig. 5a the tra- 
jectory that contains the scattering from the Friedel os- 
cillation is shorter than that of the reference trajectory, 
whereas it is longer in the configurations of Figs. 5b- 
e. The phase difference between the scattered trajectory 
and the reference trajectory is precisely compensated by 
the phase of the Friedel oscillation.^^ A similar interpre- 
tation applies to the Fock contribution, although here 
the scattering is from Friedel oscillations of the density 
matrix, not of the electron density itself. 



FIG. 5: (Color online) Configurations of classical trajectories that give the Hartree contribution to the interaction correction 
SG AA- The configurations in parts b-e contain a small-angle encounter, indicated in blue. All five configurations also contain 
a crossing of the classical trajectories. The momenta associated with the crossing are denoted pi and p2. 



E. Coulomb interaction 

For the Coulomb interaction f/c(ri,r2) = e^/|ri — r2|, 
due to the long-range nature, it is never sufficient to deal 
with the first order in perturbation theory only, and ef- 
fects of dynamical screening have to be included. Hence 
it is not sufficient to consider only diagrams with a single 
bare interaction line as in Fig. 2; instead one has to sum 
up a ladder of diagrams within the Random Phase Ap- 
proximation (RPA). This analysis is explained in Refs. 
9,26,33,40 and it can be carried over to the semiclassical 
formalism without significant modifications. 

For the purpose of including higher order interaction 
contributions, the separation into Hartree and Fock con- 
tributions is no longer meaningful. Instead, it is favor- 
able to decompose the interaction into singlet and triplet 
channels. Hereto, we consider the interaction amplitude 
of a scattering process, where two particles with initial 
momentum pi (p2) and spin a (7) interact and depart 
with final momentum Pi + q (p2 — q) and spin P (S). 
Since in the semiclassical theory of transport only paired 
trajectories are relevant, we may restrict our analysis to 
the case |q| <C Pf- The classification into singlet (j = 0) 
and triplet part {j = 1) then amounts to separating the 
interaction amplitude according to its spin structure as 



where a* represents the Pauli- matrices {i = x^y^z). To 
lowest order, such interaction processes are given by the 
Hartree and Fock scattering as depicted in Fig. 6. While 
the Fock-process is of singlet nature, the Hartree-process 



has to be split into singlet and triplet contribution. 
ut°\cD = Uciq) - liUciPl - P2))|p,|=|p.|=p., (32) 

Ut'\^) = -liUciPi - P2))|p,|=|p.|=,., (33) 

Here we anticipated, that for a diffusive system we may 
average over the directions of the momenta of the elec- 
trons. To include screening, one then sums up the RPA- 
series in each channel, 

C/«(q,t^) = U^'\ci) - [/(^■)(q)n(q,w)[/«(q,c^) (34) 

where the disorder- averaged polarization operator H is 
diagonal in spin space (note, that the disorder average 
of the polarization operator does not involve a crossing, 
and therefore has no TE-dependence). 

The discussion so far is valid for weak interaction (i.e. 
interaction parameter <C 1). For stronger interactions, 
one should also include Fermi liquid effects. Then the 
structure of the screened interaction remains the same, 
but the bare interaction is now expressed in terms of 
Fermi-liquid parameters 

Ut°\^) = Uciq) + ^Fo^ Ut'\^) = ^F^- (35) 

Let's turn back to the Altshuler-Aronov correction. 
Applying the preceding analysis, we find that Coulomb 
interaction is properly included, if we calculate the Fock- 
type diagrams as in Fig. 2, where the interaction is re- 
placed by the effective interaction 

[/(n, V2',uj) = U^^=^^ (n , r2; c^) + 3U^^=^^ (ri, r2; c^), 

(36) 
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pi,a 



P2,7 



pi,a 





P2,7 



+ 



Pi+q,/3 P2-q,5 Pi +q,/3 P2 -q,^ 



FIG. 6: (Color online) Fock (left) and Hartree (right) scat- 
tering process of two particles with initial momentum pi 
(P2) and spin a (7), and final momentum pi + q (p2 — q) 
and spin /3 (6). Since the interaction conserves spin, the 
Fock process is proportional to SapS^s and belongs to the 
singlet channel, while the Hartree process has the structure 
^aS^p-f = ^{SapS^s + cTa^^^-is) ^ud thcrcforc splits into sin- 
glet and triplet contribution. 



where the factor 3 comes from the spin summation and 
accounts for the multiplicity of the triplet channel. The 
precise relation between U (ri, Y2 ] uj) and U (q; uj) follows 
from the solution of the diffusion equation and will be 
clarified in the next section (see Eq. (48)). 



III. INTERACTION CORRECTION FOR 
ANTIDOT ARRAYS 

In this Section we apply the theory developed in the 
preceding section to quasi one- dimensional and a two- 
dimensional antidot arrays. The antidot arrays consist of 
a ballistic electron gas with randomly placed disc-shaped 
scatterers of size much larger than the Fermi wavelength. 
The classical dynamics in such an antidot array is chaotic, 
and diffusive on length scales much larger than the disc 
size a or the distance between discs. In particular, since 
the Ehrenfest time te = A~^ln(a/AF) ^ because 
of the large logarithm, and since is comparable to 
the transport time r, the diffusive dynamics applies for 
timescales down to te- 

A diffusively moving particle quickly loses its memory 
about the direction of motion, so that the classical propa- 
gators P(X, X'; t) depend on the positions r and asso- 
ciated with the phase space points X and X' only. This 
leads to significant simplifications of the general semi- 
classical expressions for the interaction correction ^Gaa- 
In order to evaluate 5Gaa in this limit, we start by ex- 
pressing the integration over momentum at fixed energy 
as an integration over the corresponding angle 0, 



dp^ = (27r/l)Vg, 



(37) 



where u is the density of states per spin. We then find 



P(X,X';i) = 



1 



P(r,r';i), 



(38) 



where P(r r';t) depends on the positions only and 
satisfies a diffusion equation. 



with diffusion coefficient D. For a rectangular sample 
of dimension L x W ^ coupled to ideal leads at x = 
and X = L and with insulating boundaries ai y = and 
y = W^ the solution of Eq. (39) reads 

P(r, r'; t) = e{t) ^ Vq(r)^q(r')e-^'''*, (40) 
q 

with the function 
V'q(r) = 



^sin(oa:)xr/^ = (41) 



The sum over q runs over = ^ with n = 1, 2, ... and 
Qy = ^ with m = 0,1,.... We also use the Fourier- 
transformed diffusion propagator. 



Do^ — iuj/h' 



(42) 



Finally, the probabilities Pin(r) and Pout(r) that a tra- 
jectory originating at position r exits the sample on the 
left or on the right, respectively, are 

Pout(r) = |, P^{v) = ^, (43) 

which may be derived from the diffusive flux 

j,(r,r^t) = -i^a,P(r,r^t), (44) 

at position r and time t, for a particle starting from at 
time t = 0. The Drude conductance (Eq. (15)) is then 
expressed as 



which gives the famihar result 



(45) 



(46) 



Let us now turn to the interaction. For the inclusion 
of screening effects we need the polarization operator, 
which, for the low frequencies at which the electron dy- 
namics is effectively diffusive, can be expressed through 
the diffusion propagator. 



n(r, r';uj) = 2v [^(r - + f P(r, ( 



(47) 



(39) 



Using Eqs. (34), (35), (36), we then find, that the effec- 
tive interaction can be written as 

C/^(ri,r2;c^) = ^ Vq(r)Vq(r')C^''(q,^) (48) 



where ?7^(q,w) = [/^'(^■=o)(q;w) + 3C/^'(^=i)(q;w) is 
given by 

!;«-l«,;„) = 153i^ (49) 
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Dq^ — iuj/h 



2v Dq^{l^ F^) -iuj/h 



(50) 



in the singlet and triplet channel, respectively. Due to 
the divergence of the bare Coulomb interaction at small 
momenta, the interaction in the singlet channel is set by 
the polarization operator solely. In the triplet channel. 



the interaction depends on the zero angular harmonic 
of F^, which is the only free parameter controlling the 
interaction strength. 



For the further calculations, it is convenient to make 
use of the diffusion equation, to write 



d 



[P(r,rVE)e'""^/^] 



With the help of additional spatial integrations over delta functions, which we then replace using Eq. (42), the 
interaction correction to the conductance then takes the form 

^Gaa = J dujduj (ujcoth-^^lml^J dr idr2U^{ri,r 2; uj) K {r uj)^ , (51) 

where the function K{ri,r2;oo) reads 

K{ri,r2;uj) = c^rc^r'P,,(r)Po,t(rO |^>c.P(r, ri; c^)p:,P(r2, c^) ^ (itP(r^ r; t)e^"*/^ 

- P(r, n ; a;)P(r2, uj)V^P{v' , r; rE)e^"^^/^ + P,P(r, n ; a;)P(r2, a;)P(r^ r; rE)e^"^^/^ 

+ P(r, vr,uj)V'^P{v2, v'- uj)P{v', r; rE)e^"^^/^| , (52) 



with the short-hand notations V^^ = {DAr + iuo/K) 
and = (DAr' + iu/h). The technical advantage 
of the structure of Eq. (52) is that each term con- 
tains the same diffusive propagators. Performing sev- 
eral partial integrations, using APin(r) = APout(r) = 
0, VPin(r) = -VPout(r) = -7;ex, as well as 
V^J^dtP{r',r;t)e''^^^^ = -P(r^ r; rE)e^^^E/^, we are 
able to simplify the expression and finally obtain 



4P)^ f 

K(ti,T2;uj) = - / dTdT'P{T,Ti;uj)P{T2,r'; 



poo 

X dxdx' / dtP{r',r-t)e 



iojt/h 



(53) 



Together, Eq. (51) and (53) represent the main result 
for the interaction correction in antidot arrays. For 
zero Ehrenfest time, the time integral of Eq. (53) equals 
P{t\t;uj) and one recovers the results for quantum impu- 
rities, obtained by standard diagrammatic perturbation 
theory (see Ref. 41, where the symbol F of this reference 
equals F/2 = —Fq /{I + Fq) and the reference misses a 
factor two for the triplet contribution). If the Ehrenfest 
time is finite, it poses a short-time threshold and only 
electrons with a travel time larger than te contribute to 
the interaction correction. We now discuss Ehrenfest- 
time dependence of SGaa in detail for a quasi- Id and a 
2d antidot array. 



A. Quasi-one dimensional antidot array 

For a quasi- Id antidot array (width W much smaller 
than length L) , we may simplify the diffusion propagator 
by taking only the diffusion mode with zero transverse 
momentum into account. After insertion of the diffusion 
propagators and the interaction into Eqs. (51) and (53), 
and using the residue technique for the cj-integration, we 
find 



2 00 00 



m=l n=l 
TeTd 



9m{nTB/TT) - n^g'^{nTj^/TT) 
T 



(54) 



where the time tt = hj^nT is the inverse temperature, 
td = jD-K^ is the diffusion time, and the function 
gm{x) is expressed as 



128^ 



k=l ^ ^ ^ ^ 



+ 



3F^ 



fc2 k'^(l + F^)+x 



(55) 
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with 



{km/(k'^ 




m2) 



if 

else. 



m odd, 



(56) 



The summation over k in Eq. (55) can be written in a 
closed form (see Eq. (Bl) in the appendix). We wih now 
discuss the dependence of 5Gxa on temperature, system 
size and Ehrenfest time several limiting cases. 



1. td ^ tt,te 

We first consider the limit td ^ oo corresponding to 
a large antidot array. In this case, the Ehrenfest- time 
dependence of the interaction correction SG aa is deter- 
mined by the ratio te/tt- For te/tt <C 1 one finds the 
result 



1 



2yr 



^0^ 



(57) 

independent of te and known from diagrammatic per- 
turbation theory. Here C{3/2) ^ 2.61238 is the Rie- 
mann zeta function (see Appendix B for details). On 
the other hand, for large Ehrenfest times or, equivalently, 
higher temperatures, te/tt ^ 1, the interaction correc- 
tion SG AA acquires an exponential dependence on tem- 
perature oc e"^^^'^^/^^ 



5Gaa 




h 7r3/2 



(58) 



independent of Fq . The crossover between these two 
limiting cases is shown in Fig. 7 for different values of 
the Fermi-Liquid interaction constant Fq . 

We emphasize the influence of the interaction constant 
in the triplet channel F^: While for small values of Fq ^ 
SGaa is always negative and monotonously decaying as 
temperature is increased, a more interesting behaviour is 
observed at large interaction strength: At small Ehren- 
fest time and Fq < the contribution from the triplet 
channel dominates and gives rise to a positive sign of the 
interaction correction. On the contrary, if the Ehrenfest 
time is large, the prefactor of the exponential behaviour 
shows no dependence on Fq to leading order in ^ and 
is therefore always negative. Hence, at sufficiently large 
interaction strengths, one observes a sign change of the 
interaction correction, as the temperature is varied. 



2. tt td, Te 

In the limit of zero temperature SG aa is a function of 
the ratio te/td only (we again refer to Appendix B for 
details). At zero Ehrenfest time, we have 



SGaa 



Li 



dxgmix), 



(59) 
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FIG. 7: (Color online) Interaction correction to the conduc- 
tance of a quasi-one dimensional antidot array in the large- 
system-size regime td ^ tt, te- 
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FIG. 8: (Color online) Interaction correction to the conduc- 
tance of a quasi-one dimensional antidot array in the low- 
temperature regime tt ^ td, te. 



with gm{x) defined in Eq. (55). For small Fq", we have 
SG AA ^ — (e^//i)(0.74 + Fq"). In this parameter range 
the singlet contribution is dominant, which leads to a 
negative sign of the interaction correction, and to a 
monotonous decay as a function of At larger Fq 
the triplet contribution competes with the singlet con- 
tribution, resulting in a sign change of the interaction 
correction for sufficiently strong interactions, starting at 
-0.5. For Te ^ td, we find an exponential depen- 



dence of SG aa on ^ 



SG AA 



h 



3F^ \ 16(57r2 - 48) td 



1 



-IE 



TE 



(60) 

The crossover between the limits te <C td and te ^ td 
is shown in Fig. 8 for several representative values of Fq . 



3. Te ^ tdjTt 

Finally, in the "classical" limit of large Ehrenfest times, 
we have 



^Gaa 



_ IE _ IE. 



(61) 



12 



where gi is defined in Eq. (55). In this parameter regime 
the interaction correction has the characteristic exponen- 
tial suppression SG aa oc e-^^^^E/^-TE/rD^ 



B. Large 2d system 

We will now consider a two-dimensional antidot ar- 
ray of dimensions L x W, where we restrict ourselves 
to the limit of large system size, td ^ te, tt- In the 
large-size limit, the relevant quantity is the conductivity 
a = GL/W. We may then express Eqs. (51) and (53) in 
momentum space. 



SCTAA 



coth 



X Im < 



(2^)2 



2TJ 



(Dq2 -iuj/hf 



(62) 



In the limit of zero Ehrenfest time, this expression sim- 
plifies to the well-known result of diagrammatic pertur- 
bation theory.^'^^ The full Ehrenfest-time dependence is 
shown in Fig. 9. For te <C tt, we have the asymptotic 
behavior 



^CFAA = T 



3 1 



ln(l + Fo-) 



^0^ 



In—, (63) 



which coincides with the well-known expression of quan- 
tum impurities, where the role of the elastic scattering 
time as a short-time cutoff is taken over by the Ehren- 
fest time. In the opposite limit te ^ tt, we obtain an 
exponential dependence on temperature. 



ScTAA T-e . 

Tin 



(64) 



As in the one-dimensional situation, at small Fq the sin- 
glet contribution dominates the interaction correction, 
while at larger Fq the triplet contribution competes, and 
a sign change of the interaction correction as a function 
of ^ is observed if F^ < -0.45. 



IV. CONCLUSION 

In this article, we considered the effect of a finite 
Ehrenfest time on the interaction correction of a con- 
ductor in which the motion of the electrons is described 
by chaotic classical dynamics. Using semiclassical theory 
of transport, we derived an expression for the interaction 
correction containing only the interaction propagator and 
coarse-grained classical propagators of the electronic mo- 
tion. We confirm the result of Ref. 19, obtained for a 
double ballistic quantum dot, that the Ehrenfest time 
enters as a short-time threshold for the interaction cor- 
rection. In other words, the minimal time it takes to 
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FIG. 9: (Color online) Interaction correction to the conduc- 
tivity cr of a two-dimensional antidot array for various values 
of the interaction strength Fq in the triplet channel. 



traverse system for trajectories responsible for the inter- 
action corrections is the Ehrenfest time. 

As a specific and experimentally relevant example 
we applied the formalism to antidot arrays, where 
the coarse-grained classical dynamics follows a diffusion 
equation. At zero Ehrenfest time, we recovered the well 
known results of the diagrammatic perturbation theory 
for a disordered metal. ^ If the Ehrenfest time is large, 
we found that the interaction correction is exponen- 
tially suppressed ex e~^^^^^ e~'^^^^^^^ . While the factor 
g-TE/TD Q^igQ present for weak localization, the sup- 
pression with temperature is specific to the interaction 
correction. Unlike the dwell time td, which governs the 
Ehrenfest-time dependence of weak localization, temper- 
ature is a variable that can be easily controlled experi- 
mentally without changing the classical dynamics, mak- 
ing the interaction correction a promising experimental 
signature of the Ehrenfest-time dependence of quantum 
transport. (We note that weak localization depends on 
temperature implicitly via its dependence on the dephas- 
ing time. However, an independent measurement of the 
dephasing time that enters into the expression for the 
Ehrenfest-time dependence of the weak localization cor- 
rection is problematic.^^) 

A particular signature of the underlying classical mo- 
tion is a sign change of the interaction corrections for 
strong enough interactions. Associated with this sign 
change is a non-monotonous temperature dependence 
of the interaction correction, in the temperature range 
T H/te- As long as only the Fock contribution is con- 
sidered, the sign of the interaction correction is negative. 
If the Hartree contribution — more precisely, the triplet 
channel of interaction — is added, there is a competition 
between Hartree and Fock- type corrections, and the sign 
of the total interaction correction at small Ehrenfest time 
may change as a function of the interaction strength. For 
large Ehrenfest time, the Fock contribution always dom- 
inates, so that the sign of the interaction correction is 
independent on the interaction strength in that limit. 

As mentioned above, the sign change of the interac- 
tion corrections for systems with small Ehrenfest times 



13 



requires a rather strong interaction in the triplet channel. 
In particular, in 2d systems the threshold was estimated 
to be Fq < —0.45, where Fq is the corresponding Fermi 
liquid constant, whereas in quasi one-dimensional sys- 
tems the condition reads Fq < —0.75. Let us focus on 2d 
systems with Coulomb interaction, for which the condi- 
tion is less restrictive. Both numerical and experimental 
results for Fq are available in the literature. In general, 
Fq is a function of the gas parameter = \f2e^ jehvp^ 
where s is the static dielectric constant and vp the Fermi 
velocity. In Ref. 42 the constant F^ was extracted from 
experimental data using the results for the interaction 
corrections to conductivity^'^'^^'^^. For systems with 
?^ 1, typical values of F^ were found to be of the 
order of —0.35. This seems consistent with numerical re- 
sults obtained in Ref. 44, were systems with moderately 
large Vg were analyzed. For the maximal value of = 5 
considered in this paper, the Fermi liquid constant de- 
creased further down to F^ = —0.5. Considerably larger 



negative values up to F^ ^ —0.7 were inferred for sys- 
tems with Ts ~ 22 in Ref. 45. We take this as evidence 
that 2d systems with sufficiently strong triplet channel 
interactions are realizable, provided the additional anti- 
dot structure may be superimposed. Relevant values of 
the gas parameter are likely in the range of ^3 — 5. 
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Appendix A: Details of the semiclassical calculation 
1. Sum rule 



Here we show how to derive the sum rule Eq. (11). We start by noticing that 



Ai 





( 9Po 












\ dr 


d^ ) 



Hence we may write 

E ^«/(P«>P«>^«) = l^'^^l'^P' E ^(*-^a)^(P'-Pa)/(P''P«' 

The determinant serves as a Jacobian for the transformation (t,p') 



t) 



( dr 


0L\ 


det 


1' 


\ dt 


dt / 



(Al) 



• (A2) 



E <5(i-T„),5(p'-p'J/(p',p„,t) 

air'— )-r;^ 



/ dr d^ 

det % % 

\ dt dt 



= S[^-Hiv,p{r',p';t)mr-v{r',p';t)] 

xf[p',p{v',p';t),t], (A3) 



where (r(r', p'; t), p(r', p'; t)) is the phase space point that a trajectory originating from (r',p') reaches after time t. 
After insertion of J dpS{p — p(r',p';t)) we finally arrive at Eq. (11). 



2. Convolution rule 



In this appendix we derive the convolution rule 

J dridr2a^(r',r2;0a^(r2,ri;^-a;)aA(ri,r;0/(ri -ra) 

1 9^ pTa rt 

= -n^JZ^^.j:j<^^-''^"'l I ^^'/(r„W-r„(0)e-(-)/^ (A4) 

where /(r) is an arbitrary function. (Here we omitted contributions from stationary configurations of trajectories, 
that cannot be connected to a single trajectory. For the calculation of the interaction correction, such contributions 
drop out upon pairing with the retarded trajectory.) 
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To this end, we first prove an auxiliary identity for the stabihty amphtudes. Hereto we consider a trajectory a that 
connects with r. Further, let ri be a point on trajectory a, and a' {a") be the part of trajectory a connecting 
with ri (ri with r). The stability amplitude can be then written as^^ 



1 


a25„(r,r') 


1 










dr± 



(A5) 



where is the Fermi velocity and r± (r^) denote displacements perpendicular to the trajectory a. The last equation 
implicitly defines r±oc{f^±^P±)- We then introduce 



dp' 



dr^gn 



la" 



such that Al 



1^-1 1, Al, = v~'^ \B~^l Al„ = v~'^ \B~^\. Then the following identity holds: 



Ba — Bot'Bot" 



aX'(ri,r') , 525„.(r,ri) 



(A6) 



(A7) 



For the proof of Eq. (A7), we note that Ba measures the change of the final coordinate of a induced by a small 
change of the initial momentum. When we consider a to be composed by a' and a small change of the initial 
momentum leads to a change of the intermediate coordinate and momentum, which results in a change of the final 
coordinate. 



dr±a 

dp±ia" 
dp±ia" 



dp^ia' 

dp'^a' 
9P±a' 



dr 



dp^la' 



P±i 



9P±a' 



dpi 



9pla' 
dp±la 



dr 



dpi 



dr 



dr^ia" 



(A8) 



The last line yields Eq. (A7). 

In a similar fashion, one verifies that 



Bol' = BaBa'i 



d^Sa"{v,v^) d^Sa{v,v') 



drl 



drl 



(A9) 



The identity we need for the derivation of the convolution rule (A4) involves the partitioning of a single trajectory 
a into three trajectories a' {v' ^ ri), a" (ri V2) and a'" {v2 v). In this case, we have 



Boi = Boc'Boc"Bo 



d^ d^\ f d^Sg. _ / d^Sg. 



9rli ^ drl^ 



dr\2 ~^ 



±2 



\dr^idri_2 



(AlO) 



To see this, one introduces the trajectory a as connection of a and a' and makes use of Eqs. (A7) and (A9). 
We now turn to the proof of the convolution rule (A4). Hereto, we define 

K(r, r'; u) = j dvidv25(vi - V2 - a)^^(r', rs; 0^''(r2, ri; ^ - c^)a^(ri, r; 0, 

where a is arbitrary, but fixed. With the abbreviation fi = ri — a we write 

if{r,r»= j drigA(r',fi;e)a^(fi,ri;^-w)gA(ri,r;C). 



(All) 



(A12) 



We then insert the semiclassical expressions for the Green functions which expresses the former equation as a sum 
over trajectories a' (from to fi), a" (from fi to ri), and a'" (from ri to r). 



15 



The integration over ri is carried out within stationary phase approximation. Here we only take into account 
stationary phase configurations, where a', a" and a!" are connected to a single trajectory. Other configurations play 
no role for the calculation of the conductance, since the three advanced trajectories are paired with a single retarded 
trajectory. Hence the convolution K may be written as a sum over trajectories a connecting with r and a sum over 
points r^^^ , for which ol first passes through f = r^^^ 
v^^^ may be parametrized as 



a and then through r^^\ Deviations Ari = (Axi, At/i) from 



Ar^i = — sin(6>)Axi + cos(6>)A?/i 
Af^i = -sin(6')Axi + cos((9)A?/i 



(A13) 



where Q (0) is the angle of the momentum of trajectory a at r^^^ (ri^^)^ and in turns Ar±i {Ar±i) represent perpen- 
dicular displacements of trajectory a at r^^^ (i^i^^). We then expand the sum of the actions of trajectories a\ a!' and 
a"' up to second order in Ari : 



<Sc/(fi,r';^) +5c//(ri,fi;<f - cj) + 5^/// (r, ri; ^) = 5c(r,r';^) - cjr + A5(Ari), 
where r is the duration of a between f^^^^ and r^^^^ and A<S(Ari) is given by 



(A14) 



A<S(Ari) =- 



+ 



(0) ;;(0)^ 



9f_Liarj_i 



AriiAf_Li, 



(0)^ 



+ 



(A15) 



where as a consequence of energy conservation only perpendicular displacements need to be considered. The integra- 
tion over Ari can then be accomplished and using Eq. (AlO) we get 



if(r,r';a;) 



27r 



(27ri/i)3/2 



«:rt:r;«^^/F|sin(^-0)| 



where the factor | sin(^ — originates from the Jacobian of the transformation (A13). (A possible phase shift from 
taking the squareroot of Eq. (AlO) is needed to restore the correct Maslov index. For our calculation however, the 
Maslov index plays no role and we drop it in our expressions.) 
On the other hand, we have 



r dt fdt'S'-'\r„{t)-r„{t')-a)= ^ 
Jo Jo /T^ Vp sm{0 



1(^-^)1 



With that, we finally obtain 

1 2-K 



y Ao.e-'^"/^ n dt f dt'S'^^\r„(t) - r„{i') - a)e'' 
^0 Jo 



,(t-t') 



{A16) 



air'— )-r;^ 



Multiplying with /(a) and integrating over a then yields Eq . (A4). 



3. Summation over classical trajectories involving a 
small-angle encounter 

The summation over classical trajectories with a small- 
angle encounter as given in Fig. 4 is performed using the 
procedure of Refs. 14,23. We here outline the main points 
of this calculation. 

The four trajectories a (from to r2), /3 (from ri to 
r), 7 (from to r), and 6 (from ri to r2) are piecewise 



paired as shown in Fig. 10. We start by noting that the 
choice of the retarded trajectories a and P fully specifies 
the advanced trajectories 7 and since the linearized 
chaotic dynamics allows for precisely one unique solution 
of a trajectory that satisfies the initial and final condi- 
tions required for the pairing shown in Fig. 10. More- 
over, the products of the stability amplitudes are equal, 
Ao,Aj3 = AjAs, so that the product of four Green func- 
tions required for the calculation of SGaa can be written 
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as 



^^^(r, n; Oa^(ri, r2, C - uj)g\v2, r'; £,)g\v' , r; i) 



air'— )-r2;^ /3:ri— )-r;^ 



(A17) 



where AS is the action difference -\- Sf3 — — Ss- 
The summation over trajectories a and /3 is restricted to 
those trajectories that undergo (at least) one smaU-angle 
encounter. 

The action difference AS has two contributions: One 
contribution from the length difference of the retarded 
trajectories a and /3 vs. the advanced trajectories 7 and 
(5, and one contribution from the different energy ^ — 
associated with the trajectory 5. The former contribu- 
tion has been calculated in Refs. 46,47 and equals SqU^; 
the latter contribution equals urs^ where rs is the dura- 
tion of the trajectory S. Note that the product SqUq is 
independent of the position of the phase space point Xe 
along the encounter. 

The trajectories a and /3 are enumerated by first pick- 
ing a phase space point Xg on the trajectory a. The 
Poincare surface of section at this point may be param- 
eterized with stable and unstable phase space coordi- 
nates, which are denoted Sq and i^e respectively. Mov- 



ing the Poincare surface of section along the trajectory, 
the unstable (stable) coordinate grows (shrinks) as e^^^, 
where A is the Lyapunov coefficient. We choose the ori- 
gin of the coordinate system such that the trajectory 
a pierces the Poincare surface of section at coordinates 
(se^Ue) = (0,0). The point Xe is part of an encounter 
formed by trajectories a and /3, if /3 passes through 
the Poincare surface of section at phase-space distance 
|se| < c and \ue\ < c, where c is a cutoff scale, be- 
low which the chaotic classical motion can be linearized. 
(One can always simultaneously rescale the coordinates 
s and li, such that the cut-off scale is the same for both 
coordinates.) The cut-off scale c enters the final results in 
the combination ln(c^//i) only, so that the precise value 
of c is unimportant, as long as c represents a scale char- 
acteristic of the classical dynamics. One verifies that 
the choice of the phase space point Xe on the trajectory 
a and of the phase space coordinates for the trajectory 
(3 also specify the two remaining trajectories 7 and S. 
Indeed, since 7 is paired with a before the encounter, 
and with (3 after the encounter, it pierces through the 
Poincare surface of section at the same unstable coordi- 
nate as P and the same stable coordinate as a. Similar 
considerations apply to the trajectory S. The summation 
over trajectories a and /3 is then written as 



/r pc roo POO pToc rTfB 

dp^dp'^dpi^^dp2,^ / d'Ke / ds^du^ / dr^ / dT/3 / dt^ / dt/3 — 

a:r-^r2;^p:ri^r;, -^-^ -^0 -^0 -^0 



X p^(r', ^ Xe; Tc, - ta)P^(Xe ^ Fs, P2; ^a) 

X p^(ri,pi ^X*;t^)p^(X* ^r,p;r^ -t^)..., (A18) 



where X* is the phase space point located at phase-space 
displacement (se, u^) from X and the dots indicate an ar- 
bitrary function of the end-point coordinates of the tra- 
jectories a and /3. The time tenc denotes the duration 
of the encounter. The factor tenc in the denominator ac- 
counts for the fact, that Xe can be chosen anywhere along 
the encounter. The ends of the encounter are determined 
by the condition that max(|s|, \u\) = c, or that one of the 
trajectories involved ends, whichever occurs first. Since 
the phase space coordinates s and u decreases/increase 
exponentially upon proceeding along the trajectories a 
and with a rate given by the Lyapunov exponent A, 
one finds that tenc is given by the expression 

tenc = min[A"Mn(c/|se|),t/3] +min[A"Mn(c/|'?/e|),ta]- 

(A19) 

The four different scenarios, depending on whether 
and ta are larger or smaller than A~^ ln(c/|5e|) and 
A~^ ln(c/|iie|), respectively, correspond to the four contri- 
butions to SG^'a^-5G^'a^ to (^Ga'a- The same situation 



also occurs in the calculation of the shot noise, see Ref. 
48,49. 



The next step in the calculation is to replace the exact 
trajectory densities by the coarse-grained ones. In or- 
der to account for correlations inside the encounter, we 
introduce phase space points X and X' at the beginning 
and end of the encounter — if the phase-space points Xi 
and X2 associated with the interaction position are not 
inside the encounter. Outside the encounter region we 
may replace the product of the exact trajectory densi- 
ties by the product of the classical propagators, whereas 
inside the encounter only a single classical propagator 
remains. After coarse graining the classical propagators 
are insensitive to small phase-space difference between 
Xe and X*, so that we may replace X* by Xe. As a 
result, the integration over and i^e and the integration 
over Xe can be performed separately. The integration 
over Xe can be performed using a convolution rule for 
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FIG. 10: (Color online) Schematic drawing of an encounter, 
formed by the trajectories a, /3, 7, and 6, with timescales and 
phase-space points as described in the main text. 



the classical propagators, 

J dXeP(X2,Xe;t2)P(Xe,Xi;ti) = P(X2, Xi ; ti + ^2)- 

(A20) 

The contribution SG^^^ is then expressed as Eq. (22) 
with replaced by 

i^2a(Xi,X2;w) 

dXdX'Pi„(X)P(X, Xi; u;)P(X2, X'; w)Pout(X') 



I' 



X [ dsduP{X\X]tenc)e'' 

27Th J_^ tenc 



^iujtezic/hl 



, (A21) 



where P(X, X';cj) is the Fourier transform of 
P(X,X';t), see Eq. (29). The encounter time for 
this contribution is given by tenc = ln(c^/|sii|). In 
order to perform the integration over the phase-space 
coordinate s and u we make use of the integral identity 



1 

2^ 



oisu/h 



dsdu 



d 



/(tenc) = ^/(te), (A22) 

^enc 



which holds for an arbitrary function f(t) which is a slow 
function of its argument on the time scale A~^. In this 
equation, the Ehrenfest time is defined as 



TE = A-Mn(cV^). 



(A23) 



One notes that this definition is consistent with Eq. (1) of 
the main text, since is a classical action characteristic 
of the classical motion, which can also be expressed as 

= pf<^, where a is a length scale characteristic of the 
classical motion. The equivalence to Eq. (1) then follows 
since = 27rh/XF. With the equality (A22) the result 
(25) of the main text follows immediately. 

To prove Eq. (A22), one makes use of the identity 



2 r 

Jo 



smx 
dx = 6>(-t), 



(A24) 



where the Heaviside step function 0{t) is broadened on 
the scale A~^. Taking derivatives on each side, we obtain 



2A2 



^sin(e-^*) = 5it), 

_d_ 
~ dt 



cos (e-^*) 



(A25) 
5{t). (A26) 



These equations can be apphed to the left hand side of 
Eq. (A22) after writing the integration in terms of pos- 
itive s and u and after a variable change that uses ienc 
and w = s/c as the integration variables, 



^g — A(tenc— Te) 



— — / dsdu f{ts^tu) 
2 

^ Jo Je-^*en 

^enc 

= — / dtence-^(^----)cOs(6 

^ Jo 



-A(tenc-rE) 



^/(te), 



)f{tenc) 

(A27) 



where we used Eq. (A26) in the last line. 

F 26 

For the contribution SGj^j^ , the encounter is bounded 
to the right by the phase space point X2 of the inter- 
action. In this case, the encounter time is given by 
^enc = A~^ ln(c/|5|) + to;, whcrc ta cau take values be- 
tween zero and A~^ ln(c/|2i|). The expression for SG^^ 
is again of the form (22), with replaced by 

^25 



K^^(Xi,X2;c^) 



dXPin(X)P(X,Xi;a;)Pout(X2) 



< J dsdu J 



27Th 



X-^ ln(c/|n|) 



dta 



xP(X2,X;tenc)e 



ia;tenc/^_ 



(A28) 



The integration over s and u in Eq. (A28) is done with 
the help of the identity 



1 



nX ^ln{c/\u\)u ^isu/h 

. dsdu / dta 

27rh Jo ^enc 



/(tenc) = /(te), 



(A29) 

which is proven by first performing the integrations over 
s and and then using the identity (A25). The final 
result is Eq. (26) of the main text. The derivation of Eq. 
(27) of the main text is similar. 

Finally, for the fourth and last contribution ^Gaa^' 
encounter is bounded by both phase space points Xi and 
X2 of the interaction vertices. The encounter time is here 
given by tenc = ta + where ta and vary between 
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and A ^ ln{c/\u\) and between and A ^ ln(c/|5|), re- 
spectively. The expression for SG^l^ takes the form (22), 
with replaced by 

= - I dXPin(Xi)Pout(X2)^ I dsdu 

n\-Hn{c/\s\) nX-Hn{c/\u\) 
X dt/3 

Jo Jo 



dta 



te 



(A30) 



The integrations over s and u can be performed with the 
help of Eq. (A24) and yield Eq. (28) of the main text. 



Appendix B: Details of the discussion 

In this appendix we add some technical details to the 
discussion of Sec. III. The function gm{x) of Eq. (55) can 
be cast in the following closed analytic expression 



where terms with an additional factor ^/ttJty) cx \ cor- 
respond to finite size corrections and will be neglected. 

After replacing the summation over m in Eq. (54) by 
an integration, we obtain 



with 

h{x) = f 
Jo 



4e2 
irh 



n=l 



-) 

TT J 



5/2 



(B3) 



oo 2 

dz 



1 



z(\^F^)^x_ 



(B4) 

Apart from the prefactor \Jtt Ity)^ the interaction cor- 
rection ^Gaa is a function of the ratio t^/tt only. The 
limiting behavior for small and large ratios te / tt is given 
in the main text. 

For the case tt ^ te , td , we replace the summation 
over n in Eq. (54) by an integration and find 



9m{x) 



32m^a tanh^ ( a 



^2 (^^2 _^ TT^y/x {m? + x) {m? + ax) 



7r2 {m? + xy ((1 + Fg)m? + x) 



966^/2^^2 tanh" 



TT^V^ {m? + x) {hm? + x) 



6=1 



a^O 



(Bl) 



E 



^2 IE 
TD 



dx{l + ^x)gm{x)e 



TE . 



' TD 



(B5) 



This is a function of the ratio te / td only. 

For the case of a large two-dimensional antidot array 
we use the residue technique to perform the uj integration 
of Eq. (62), and obtain 



where s = ±1 for m even (odd). In the regime td ^ 
tt,te, we may write this as 



gm{nTB/rT) 



1 



"^2^2 (^212: _^^)2 



1 + 



m^^{l^F^)^n 
(l + ^(V^))' 



(B2) 



SGaa 



h L TT ^ ^ 

n=l ^ 



/2 (n^) - (n^)] , (B6) 



with 

/2 



{x) = r 

Jo 



dz 



{z + xy 



1 + 



SF^z 



z{l + F^) + x 



(B7) 
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